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We discuss Cosmic Microwave Background constraints on tlie causal set theory of quantum grav- 
ity, which has made testable predictions about the nature of dark energy. We flesh out previously 
discussed heuristic constraints by showing how the power spectrum of causal set dark energy fluc- 
tuations can be found from the overlap volumes of past light cones of points in the universe. Using 
a modified Boltzmann code we put constraints on the single parameter of the theory that are some- 
what stronger than previous ones. We conclude that causal set theory cannot explain late-time 
acceleration without radical alterations to General Relativity. 



I. INTRODUCTION 

Causal sets or causets are an approach to quantum 
gravity in which spacetime is discretized with causal- 
ity as the fundamental relation between the discrete ele- 
ments. The requirement that the as-yet unknown theory 
of causet dynamics should give rise to a smooth classical 
spacetime imposes strong restrictions on the character of 
the theory. In particular, the requirement that physics 
should appear to be Lorentz-covariant implies a Poisson 
relation between the classical volume of a region of space- 
time and the number of discrete elements in it. The mean 
number of elements per Planck volume, a^, in the Pois- 
son process is the only free parameter in the theory on 
macroscopic scales. 

It was noted in Sorkin Q that the fluctuations in 
the Poisson process naturally give an appropriately sized 
dark energy if the volume of the universe is taken as 
conjugate to A. In Ahmed ct al. it was shown that 
the same principle predicts a randomly fluctuating f^A of 
roughly constant amplitude a throughout the history of 
the universe, with the name everpresent A given to the 
phenomenon. 

A fundamental objection to this explanation of dark 
energy was raised in Barrow which noted that since 
everpresent A is sourced in past light cones (PLCs) of 
points, regions of the last-scattering surface separated 
by more than a causal horizon should have uncorrelated 
A, and should produce o{a) CMB anisotropics. A limit 
of ~ 10~^ was therefore placed on a, too small for it 
to explain late-time acceleration. In this paper we flesh 
out the objection in Barrow with a numerical analy- 
sis based on the volumes of overlapping PLCs and place 
more constraints on a. 

To emphasize that we are not considering a cosmo- 
logical constant we will break from the notation used in 
Ahmed et al. [l| and use A to represent the causet dark 
energy, rather than A. 



II. THEORY 

We make use of a formulation for Everpresent A intro- 
duced in Ahmed et al. [J]: we consider the A to be the 



total Hamiltonian action S per unit spacetime volume V ^ 
with each causet element having randomly signed action 



\^S/V 



(1) 



The volume V is of the past light cone (PLC) of any 
point, and in a discrete spacetime is proportional to the 
number of elements N in it. In conformal time r, V is 
given by 

V{t) = a^N (2) 
= m\r-T'fa\T'W. (3) 

Since the action for each element is randomly chosen 
as ±a^, the total comes from binomial statistics: 



(5) = 0, 

The variance of A is therefore given by 

1 

^ ' N V 
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Since V~ 2 oc oc this formulation provides a density 
which oscillates randomly with amplitude apcrit at every 
epoch — an everpresent A. 

We can also use this formulation to determine the two- 
point correlations of A. Since A is sourced by the action 
in a PLC, we split the PLCs of two points into shared 
and disjoint regions. We use the notation A^ = \{xa) 
and consider the volume splitting that is illustrated in 
figure [TJ 
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FIG. 1: A schematic of the PLCs of two points. The Big Bang 
is at the bottom of the diagram; the dimensions are heavily 
distorted. The volumes of the light cones are Vi = Va + Vc 
and V2 = Vb + Vc. 



radiation anisotropy (see Seljak and Zaldarriaga Q). We 
use: 

STik,T)=g^, (10) 

where g = ke~'^ is the visibility function, n is the optical 
depth to time r, and -0 is the conformal gauge pertur- 
bation variable, g is computed in a standard way from 
the cosmic reionization history. We have neglected the 
effect of the metric fluctuation on the plasma, since the 
latter does not have time to respond to the former, which 
changes over a Hubble time. We have also neglected the 
Integrated Sachs- Wolfe effect. We can find an expression 
for ij} from the Einstein equations, but it will always be 
in terms of unknown stress-energy components. 

We have no stress-energy tensor for A that we could 
use to fully compute ip — we only have one component. 
We will discuss two assumptions that can be made about 
A dynamics to make progress. Implementations of each 
arc discussed in section IIIII 



The second line follows because the volumes cancel, the 
third because disjoint volumes have zero correlation, and 
the fourth by construction of S. 

We specialize to the case where the points' PLC vol- 
umes Vi and V2 are equal — the points are in a sta- 
tistically homogenous universe and have the same time 
co-ordinate. At any point in the history of the universe 
the instantaneous three- volume of the overlapping region 
is the overlap between the two points' (spherical) horizon 
volumes. The 4- volume can then be found by integrating 
this quantity. The overlap volume of two spheres of equal 
radii r and separation d is : 

V}'^^^iir + d){2r^df (7) 

if d < 2r, and zero otherwise. 

In our case for a point at time r, the horizon of a 
point a time r' earlier grown forward to t is (r — r'). For 
a comoving separation x the proper separation is a(r)x 
and the comoving horizon is r, so converting to a proper 
3- volume: 

Vj'Hr, r\ x) - a3(r)^ (4(r - r') + x) (2(r - r') - xf 

(8) 

if X < 2(r — r'), and zero otherwise. 

To compute the 4- volume of the PLC at t we integrate 
this quantity with the time element adr: 

Vo{r,x)^^f^ '''\\r){A{T-r')+x) 

.(2(T-r')-x)dr'. (9) 

We use the Vi and V2 from equa tion [3l and get the 
covariance by multiplying by (X^) — SnGap. 

To obtain CMB spectra from the A overlap we require 
a source term that converts from metric fluctuation to 



A. Zero Velocity Assumption 

If we substitute for the spatial conformal perturbation 
(f> in the perturbed Einstein equations Q we can obtain 
ip. Letting: 

D EE AnGa^ST^ = , 
T = i47rGa2(5Fr° = unkown , (11) 

then: 

k^{^ + H^) = T, 
k^(j) + 3H{^ + HiP) = D . (12) 

One approach is to neglect T altogether, and have the 
CMB sourced by A alone. This might be expected to 
give an answer of the correct order, since velocity and 
density perturbations are typically the same magnitude, 
and unlikely to cancel exactly. In that case: 

= -^k-^a^X. (13) 

where in the last line we have set A = because it is 
formally undefined. 

In this approximation we therefore take ipk = 
-a^Xk/k\ 

B. Zero Shear Assumption 

Another reasonable assumption that we could make 
would be that like many other forms of dark energy, the 
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Everpresent A has zero anisotropic metric shear. In that 
case i/j = (p and we can find tjj in terms of ST^ = A/SttG: 



which has the solution: 
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(14) 



(15) 



(16) 



We have set the constant of integration to zero so that 
there is no perturbation at zero time. The H part of the 
integrating factor can be spht off and integrated to log a 
so that: 



1 



I M Qk{r.r')^^\,{T')dT' , (17) 
2a(T) 7o H{t') 



where: 



Q,(r,T') =exp (^-h^ £ H-\T")dT" 



(18) 



Here Q is playing the role of a kernel, defining the length 
of time over which A contributes to ip. It peaks to unity 
at r' = T and falls off over about a Hubble time, or longer 
for longer wavelengths. We can compute Q simply from 
the background dynamics of the system. 

Equation [T7] relates the dynamic A and ip. For our pur- 
poses we need the relation between the root-mean-square 
quantities. This becomes somewhat more computation- 
ally involved: 



H{n) H{t2) 
(AUn)Afc(r2)) 
dridr2 . 



(19) 



We use the square-root of this quantity in the CMB 
sources. The quantity (A^(Ti)Afc(T2)) is an unequal time 
correlator similar to ones used in studies of textures and 
strings . We need to calculate it for each pair of times 
Ti, T2 where Q is not negligible, for each wavenumber k. 

As with the equal time version (to which it should 
reduce when ti = T2), we use equation [H] to compute this 
correlator. This time, though, Vi 7^ V2. To calculate V12 
we once again consider overlapping spheres, this time of 
unequal radii. The overlap volume of two spheres of radii 
ri and r2 with separation d is Q: 

K''^ = ir,+r2-df{d' + 2d{r,+r2)-3{rl+rl) + 6r,r2). 

(20) 



Once again we substitute d{T') = a{T')x, ri = 
a{T'){Ti — t') and 7^2 = a{T')(T2 — r'), and integrate: 

<.(ti+T2-x)/2 

Vo{r,,r2) = - (ri+r2-2r-x')' 

[x' + 2(ti +T2- 2t')x 

+6(ri-T')(T2-r')] 

(21) 

We get Vi and V2 from equation [H and go fr om t he 
correlation to the covariance by multiplying by (A^) ~ 
SirGap. 



III. IMPLEMENTATION 

We modified the CMB Boltzmann code CMBEASY 
3 to include the perturbations described. The code, 
which is restricted to flat universes, first computes the 
background evolution of the cosmos, which we leave un- 
changed, and then evolves perturbations (which we skip 
for the causet A) , and finally computes CMB sources and 
calculates their Ce. We compute the A contribution at 
the same values of r and comoving wavelength k as the 
code uses in its standard mode, since these are the pe- 
riods when the visibility g is significant enough for the 
perturbations to make a contribution to the CMB. 

We obtain quantities in harmonic fc-space using a Fast 
Fourier Transform routine [1] from the real space quanti- 
ties that we calculate using the volume overlap functions. 

The expression for the PLC volume in equation [3] can 
be split up into: 



Vir) 
where: 



An 



[r^Ioir) - 3r2/i(r) + 3t/2(t) - h{T)] 



In - / r'"a4(T')dr' , 
/o 



di, 
d^ 



a^r) 



(22) 

(23) 
(24) 



And thus: 



V = 



3t^Iq + T^io - 6r/i - 3r^/i 



27 



-3/2 + 3t/2 - /f 



= f (3t^/o + T^a^ - 6r/i - 3T'^a^ 

+3I2 + Sr^a"^ - T^a"^) 
= 47r (r^/o - 2t/i + 3/2) . 



(25) 



We can split the overlapping volumes of the PLCs in 
equation [5] in a similar way. The integrals range from 
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to the time when the horizon spheres are smaller than 
the comoving separation of the points. We therefore use 
the quantities: 

Jn{T,x) = r T"'a\T')dT- r T'"a\T')dT 
Jo 

= / r'"a4(r')dr, (26) 

provided the integral range is positive. We can use the 
fimctions /„ defined in equation 1241 to express J„: 

J„(T,x)=/„(r-x/2) (27) 

if T > x/2, and zero otherwise. 

Using these quantities the overlap volume is: 

+ 12(4t-x)J2 
- 24(2t2-tx)Ji 

+ (16t3-12t\ + x')Jo] • (28) 

We use this formula, with the volume formula, to com- 
pute the Afc covariancc for the CMB sources. 

To compute the unequal-time overlaps in equation!^ 
we define yet another set of retarded volumes. In this 
case the latest time r at which the horizons overlap is a 
function of both ti and T2 : 

K = ^[ - 16if3 

+ 12{xix~s) + d^)K2/x 

- 12{{x- s){sx~d^))Ki/x 

+ ix- s fix^ ~ + 2sx)Ko /x] , (29) 

where 

i^n(x,Ti,r2) ^ /„((ti +T2 - x)/2) (30) 
if Ti + T2 > Xj s-iid zero otherwise, and 

S = Ti + T2 , 

= Tl - T2 . (31) 

Since computing the overlap volumes was somewhat 
computationally intensive, and was repeated for each fc, 
we wrote code to perform it once for all the relevant k 
and pairs ti and T2 during the peak of recombination, 
and save the data. We then looked up the result in the 
main code where necessary. 

We use two other devices to speed the code: we com- 
pute the correlators only when the kernels Qk are not 
small (that is, only for ti,T2 in the recent past of r), and 
we compute the overall sources only when the visibility 
is greater than 1% of its maximum. 



IV. RESULTS AND COMMENTARY 

Our results are in the form of standard CMB power 
spectra, for temperature only since we do not consider 
any polarization sources in the calculation. In each case 
we obtain a limit on the fundamental parameter a of the 
causal set theory. 

The result from the zero velocity assumption, for a — 
10~^, is shown in figure [H It has two expected features: 
power-law increase with scale, which occurs as the the 
separation between points on the LSS increases and their 
shared volume rapidly decreases; and a kink (. k, 200, 
which corresponds to the causal horizon on the last scat- 
tering surface, and is thus where volumes become com- 
pletely disconnected. 
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FIG. 2: The power spectra obtained from a standard CDM 
parameter combination described above (red, dashed) and 
from the causal set fluctuations (black), in the case of zero 
velocity, for a = 10~®. 

Figure [3] shows the result of the zero shear computa- 
tion, for a = lO^*". The kink at the causal horizon is 
retained, but the slope at large scales is reduced because 
at large scales the terms with the fc^ factor in equationfMl 
(which were implicitly all we included in the zero velocity 
assumption) lose dominance. 

V. CONCLUSIONS 

Under both our assumptions we can place strong lim- 
its on the value of the fundamental causet parameter a, 
which controls the amplitude of all the causet fiuctua- 
tions and thus (quadratically) controls the power spectra 
scales. In both cases, and probably in all simple models, 
the causet C( are largest at low £. This makes com- 
paring the predictions to data straightforward; we can 
simply see at what a the low-£ causet Ci is equal to the 
measured value. A higher a would lead to higher Cg and 
is therefore ruled out. 

In both cases, as predicted by Barrow, the large scale 
CMB values limit a to well below 1. The exact con- 
straint depends on the detail of the conversion from the 
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FIG. 3: The power spectra obtained from the standard CMB 
parameter combination (red, dashed) and from the causal set 
fluctuations (black, solid), in the case of zero shear, for a = 
10"®. 



A perturbation to '0. In the case of the zero velocity as- 
sumption described above the limit is is approximately 
a < 10~^. In the case of the zero shear assumption the 
limit is approximately a < 10"''. In both these cases, 
again as Barrow notes, the given a is too small to gen- 



erate late time acceleration. We would require a to be 
at least o(10~^) to achieve that. We therefore concur 
with Barrow and find no loophole in the mechanics of 
causal set structure formation that would invalidate his 
conclusion that Everpresent A is not dark energy. 

From a positive perspective, we can think of these re- 
sults as providing direction to the theory. A successful 
version of it must resolve the problem of excess structure 
formation, but retain acceleration at late times. This 
might come about by finding a very unusual stress-energy 
tensor, by suppressing A at early times, or by modifying 
the Einstein equations more dramatically. The last seems 
the most promising possibility. 
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